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Nuclear attenuation of the two-hadron system is considered in the string model. The two-scale 
model and its improved version with two different choices of constituent formation time and sets of 
parameters obtained earlier for the single hadron attenuation, are used to describe available expe- 
rimental data for the z-dependence of subleading hadron, whereas satisfactory agreement with the 
experimental data has been observed. A model prediction for ^-dependence of the nuclear attenua- 
tion of the two-hadron system is also presented. 
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I. INTRODUCTION 

Nuclear Attenuation (NA) of high energy hadrons 
is a well known tool for investigation of an early stage of 
the hadronization process, which can not be described 
in the framework of the existing theory of strong in- 
teractions (perturbative QCD), because of major role of 
"soft" interactions. Nevertheless, there are many phe- 
nomenological models, which describe, rather quali- 
tatively, existing experimental data for single hadron 
NA 1 1 ]— 1 10] . Also some predictions for the attenuation 
of multi-hadron systems leptoproduced in nuclear mat- 
ter in the framework of the string model |11]-|12] were 
done. It was argued that measurements of NA of multi- 
hadron systems can remove some ambiguities in deter- 
mination of the parameters describing strongly interact- 
ing systems at the early stage of particle production: for- 
mation time of hadrons and cross-section of the inter- 
mediate state interaction inside the nucleus. Recently 
for the first time, data 1 13] on the two-hadron system 
NA ratio measured as a function of relative energy of 
the subleading hadron has been published. Therefore, 
in this work we attempt to describe these data based 
on the Two-Scale Model (TSM) |4] and Improved Two- 
Scale Model (ITSM) 1 10]. We present also predictions for 
the ^-dependence of two-hadron system NA within the 
same model. Also, possible mutual screening of had- 
rons in string and its experimental verification is dis- 
cussed. 



process can be expressed in form: 

Rfj = 2da A {v,Q 2 ,z 1 ,z 2 )/Ada D {v,Q 2 ,z 1 ,z 2 ) , (2) 

where da a and dan are the cross-sections for the reac- 
tion (1) on nuclear and deuterium targets, respectively, 
v and Q 2 denote the energy of the virtual photon and 
square of its four-momentum. One can picture the reac- 
tion (1) as shown in Fig.Q] The interaction of the lepton 
with nucleon occurs at the point (b, x), where the inter- 
mediate state q begins to propagate (6 and x are impact 
parameter and longitudinal coordinate of DIS point). At 
some points the string breaks, and as a result the first 
constituents of hadrons h\ and h 2 are produced at the 
points (6, xi) and (&, x 2 ). Also, at the points (6, x y i) and 
(6, x y2 ) the yo-yo of the hadrons hi and h 2 is formed (by 
a "yo-yo" system we assume, that the colorless system 
with valence content and quantum numbers of the final 
hadron is formed, but without its "sea" partons). 

In the string model there are simple connections be- 
tween the points x y \ — X\ = z\L and x y2 — x 2 = z 2 L, 
where L is the full hadronization length, L = v / ' n, n is 
the string tension (string constant). The NA ratio can be 
presented in the following form: 

1 /> />oo />oo />oo 

Rm ~ o / ^ / dx dx ^ / dx 2P(b, X) (3) 

J J— oo J x J xi 

[D(zi,z 2 , xi — x,x 2 ~ x)W (h 1 , h 2 ; b, x, x\, x 2 )+ 



II. THEORETICAL FRAMEWORK 

In articles ] 1 1 ]— 1 12] the process of leptoproduction 
of two-hadron system on a nucleus with atomic mass 
number A was theoretically considered for the first time: 

k + A->l f + h x +h 2 + X , (1) 

where the hadrons hi and h 2 carry fractions z\ and 
z 2 of the total available energy The NA ratio for that 



D(z 2l zi, x\ — x, x 2 - x)Wo(h 2 , h%;b, x, xi 7 x 2 )] , 

where D(z%, z 2 , 1%, l 2 ) (with 1% < l 2 ) is the distribution of 
the formation lengths l\ and l 2 of the hadrons and p(b, x) 
is the nuclear density function normalized to unity. Wo 
is the probability that neither the hadrons hi, h 2 nor 
the intermediate state leading to their production (ini- 
tial and open strings) interact inelastically in the nuclear 
matter: 

W (hx,h 2 ;b, x,xi,x 2 ) = (1 - Qi - Si- (4) 
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(H x + Q 2 + S 2 + H 2 - H X (Q 2 + S 2 + H 2 ))f A -V , 

where Qi and Q 2 are the probabilities for the initial 
string to be absorbed in the nucleus within the inter- 
vals (x, xi) and (x y i,x 2 ), respectively. Si (i = 1,2) is 
the probability for the open string containing the first 
constituent parton for hi to be absorbed in the nucleus 
within an interval (xj, x y i), and Hi (i = 1, 2) is the prob- 
ability for the hi to interact inelastically in the nuclear 
matter, starting from the point x y i. The probabilities 
Qi, Q 2l Si, S 2 , Hi, H 2 can be calculated using the gen- 
eral formulae: 

P{x,nin,X max ) = <T P p(b,x)dx , (5) 

where the subscript P denotes the particle (initial string, 
open string or hadron), up its inelastic cross section on 
nucleon target, and x TO j n and x max are the end points of 
its path in the x direction, as is shown in Fig.Q] 

III. EXPERIMENTAL SITUATION 

Recently the HERMES Collaboration has obtained, for 
the first time, data on double hadron attenuation 1131 . 
The following double ratio for leading and subleading 
hadrons has been measured: 

R$(z 2 ) = (d 2 N(zi, z 2 )/dN{zi)) A / (6) 

(d 2 N(zi,z 2 )/dN(zi)) D , 

where Zi = Ei/ v, Ei is the energy of i-th hadron, A and 
D marks that interaction takes place on nuclear and deu- 
terium targets, d 2 N(zi,z 2 ) is the number of events with 
at least two hadrons (leading and subleading hadrons 
in one event). According to the experimental conditions 
leading hadron must have zi > 0.5, and dN(zi) is the 
number of events with at least one hadron with zi > 0.5. 
In this experiment double ratio Bjfy(z 2 ) was considered 
as a function of z su bi ea ding = z 2 . Two sets of experimen- 
tal data are presented: 

1. the leading-subleading combinations ++, , +0, 

0+, -0,0-, 00 only; 

2. all hadron pairs except those with invariant mass 
near p°. 

In our work, we will use the first set, because, in our 
opinion, it contains less contribution from hadrons pro- 
duced in diffractive and diffractive dissociation pro- 
cesses. The relevant region for the investigation of the 
two-hadron system NA ratio as a function of z 2 is the 
one between z 2 = 0.1 and z 2 = 0.4. At z 2 < 0.1 the contri- 
butions from the slow hadrons coming from the target 
fragmentation become large. At z 2 > 0.4 it is difficult 
to distiguish leading and subleading hadron because in 
this region zi ~ z 2 w 0.5. 



Let us also consider the following question - to what 
extent are these data free from contribution of diffrac- 
tive p° mesons? For events dN(zi) corresponding cor- 
rections were made by experimentalists. Here we dis- 
cuss situation with d 2 N(zi, z 2 ) only. At a first glance the 
diffractive p° mesons do not give contribution in these 
events because of special conditions at choice of pairs of 
hadrons. They ensure that both hadrons in pair can not 
be produced from decay of a one p° meson, but do not 
forbid that one of the hadrons in the event was produced 
as a result of breaking of a diffractive p° -meson. Then, 
another hadron could be produced after the p° -meson 
final state interaction. 



IV. RESULTS AND DISCUSSION 

We have performed calculations for R 2l }(z 2 ) in the 
framework of TSM and its improved version. The ra- 
tio is a function of three variables zi , z 2 and v, but 
for comparison with the available experimental data we 
present R 2, j as a function of z 2 only, after integration 
over other variables according to the experimental con- 
ditions (0.5 < zi < 1 - z 2 , 7 < v < 23.5GeV,p h > 
lAGeV/c). In this paper we also present model predic- 
tions for R 2 ^ as a function of v, while performing inte- 
grations over zt, z% in the corresponding kinematic re- 
gions. They are chosen to be close to the conditions of 
the available experimental data (0.5 < zi < 0.9,0.1 < 
z 2 < 1 — zi). The upper limit of Zi < 0.9 relates to 
the case of two hadrons in final state, as well as for the 
single hadron this limit has to be equal to unity. All the- 
oretical curves presented here were obtained with an as- 
sumption that the final hadrons are pions. We used two 
expressions for Constituent Formation Time (CFT). As 
a first expression for CFT the equation (4.1) from fl2Ti 
(marked as CFT1) was used. Second expression was 
taken in a form according to Ref . 1 14], which means, that 
as CFT for the first produced constituent quark of the 
first hadron we use r c i = (1 — z± — z%)L, and for the 
first produced constituent quark of the second hadron 
- t c2 = (1 — z%)L (marked as CFT2). Three sets of 
parameters (including nuclear density functions) corre- 
sponding to the minimum values of \ 2 from Tables 1, 
2, 3 Ref. 1 10] were used for calculations. The value of 
string tension was fixed at k = 1 GeV/fm. Results for 
the double ratio R 2 ^} as a function of z 2 are presented on 
Fig. [5] On panels a), c), e) three options of the theoretical 
curves are shown: solid curves correspond to the TSM 
with CFT1; dashed curves correspond to the ITSM with 
CFT2; dotted curves correspond to the ITSM with CFT1. 
According to the ideology of the string model, the trans- 
verse size of the string is much less than longitudinal 
one. It means that the hadrons produced from the string 
have close impact parameters, and could partly screen 
one another, which in turn must lead to the weakness 
of NA (partial attenuation). To study this effect and to 
compare with the basic supposition, that two hadrons 



3 




FIG. 1: Leptoproduction of two-hadron system on nuclear target. Details see in text. 



attenuate independently (full attenuation), we consider 
partial attenuation in it's extreme case, when two had- 
rons fully screen one another, and as a result the two- 
hadron system attenuates as a single hadron. The re- 
sults of calculations within this conditions are shown in 
panels b), d), f). The two-hadron system will attenu- 
ate as a single hadron also when the two final hadrons 
appear as a result of breaking of one of the resonance. 
For instance, combinations +0(— 0) and + (0-) can 
be obtained as products of decay of the p + (p~) mesons 
produced via fragmentation mechanism in the nucleus, 
when the decay occured outside of the nucleus. Com- 
parison with the experimental data for z 2 -dependence 
shows that difference between versions is smaller than 
the experimental errors, consequently, different versions 
of model can not be distinguished by means of compar- 
ison with these data. Calculations with full and partial 
NA give close results. Theoretical curves quite satisfac- 
tory describe data for nitrogen. In case of krypton and 
xenon targets the situation is more ambiguous. While 
three middle points are described satisfactory, two ex- 
treme points corresponding to lower and higher values 
of z 2 are in much worse agreement. In our opinion, that 
possible reason could be that our model does not con- 
tain necessary ingredients for quantitative description 
of these points. 

Let us briefly discuss mechanisms which, we believe, 
give considerable contribution in the extreme points dis- 
cussed above, but are not included in present model. 



First mechanism which can lead to the increasing of the 
NA ratio at lower value of z 2 is that the part of sublead- 
ing hadrons in nucleus are protons, which are produced 
in abundance at small z, and in this region they have 
value of NA ratio larger than unity. We will try to esti- 
mate the contribution of protons in Appendix 1. Second 
mechanism which can lead to the increasing of the NA 
ratio at lower z 2 is the rescattering of produced hadrons 
in the nucleus, in result hadrons spend part of theirs 
energy for production of slow hadrons. Consequently, 
more slow hadrons arise in nucleus, than in deuterium, 
and, despite absorption, the multiplicity ratio in this re- 
gion can become close to and even larger than unity. 
Our model does not take into account the final state 
interactions of the produced hadrons, consequently, at 
present we can not calculate or estimate contribution of 
this mechanism. Concerning second extreme point at 
higher value of z 2 , which is equal z 2 = 0.44, we suppose 
that the double hadron attenuation ratio in this point is 
on the order of unity because there are two more addi- 
tional mechanisms, which are not included in present 
model. The first one had to do with the pairs of pions 
appearing in a result of breaking of coherently produced 
diffractive w-mesons, for which the coherent cross sec- 
tion depends proportional to the atomic mass number 
as A 2 . As a result, the NA ratio for heavy nuclei raises. 
It is very difficult to estimate the contribution of this 
mechanism without implementation of additional free 
parameters. Second mechanism is connected with the 
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FIG. 2: Double ratio i?f J as a function 0/22. The explanations for the points and curves are written in the text. 



smallness of integration region over z\, which in case 
of z 2 = 0.44 is equal to 0.06 l . The NA ratio for two- 
hadron system is proportional to the ratio of integrals 
over z\ on nucleus and deuterium. Taking into account 
Fermi motion or nucleon-nucleon correlations can lead 
to the extension of the integration region in nucleus and 
in a result to the increasing of NA ratio (see details in 
Appendix 2). The model gives close results for the two- 
hadron system NA ratio for the krypton and xenon tar- 
gets. 

Fig. |3] shows the prediction of model for in- 
dependence of double ratio R 2f } for nitrogen, krypton 
and xenon targets. On the panels a), c), e) three varieties 
of the theoretical curves are shown: solid curves corre- 
spond to the TSM with CFT1 ; dashed curves correspond 
to the ITSM with CFT2; dotted curves correspond to the 
ITSM with CFT1. On panels b), d), f) are shown the 



same curves calculated with additional condition that 
only first produced hadron attenuates (partial attenua- 
tion). Easy to see, that curves corresponding to the full 
and partial attenuation have different behavior at low 
values of v. In Fig. |3| it is shown, as example, case of 
krypton only. Curves marked as in Fig. [3] Lower curves 
correspond to the case of full attenuation and upper 
curves correspond to the case that only first produced 
hadron attenuates (partial attenuation). The measure- 
ment of NA ratio in the region of v from 3GeV to WGeV 
is allowed to verify a supposition about possible mutual 
screening of hadrons in string. We think that such exper- 
iment can be useful for comparison with the results ob- 
tained at RHIC by STAR Collaboration 0, which state 
that two hadrons from one jet absorbed more weakly 
than two hadrons from away-side jets. 



V. CONCLUSIONS 



this value of 0.06 is in fact related to the case of only two hadrons in 
final state, more hadrons in final state lead to the decreasing of the 
integration range 



• String model till Il2ll gives natural and simple 
mechanism for description of the two-hadron sys- 
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FIG. 3: Double ratio as a function off. The explanations for the curves are written in the text. 



tern NA, which allows to describe the available ex- 
perimental data for z-dependence of subleading 
hadron on a satisfactory level, using the sets of 
parameters obtained in Ref. 1 10] for single hadron 
NA. 

• Comparison with the experimental data for z-r 
dependence show that difference between ver- 
sions of model is smaller than experimental er- 
rors, consequently, they can not be distinguished 
by means of comparison with these data. 

• Double ratio considered as a function of partial en- 
ergy of subleading hadron z-i has weak sensitivity 
to the mutual screening of hadrons 

• It is of certain interest to also study other aspects 
of the two-hadron system production in nuclear 
medium. In particular we propose to measure the 
^-dependence of NA, because, as we have shown, 
it is more sensitive to the mutual screening of had- 
rons than the ^-dependence. Investigation of the 



^-dependence in the region of v from 3GeV to 
WGeV will allow better understanding of ques- 
tions connected with possible mutual screening of 
hadrons in the string. Corresponding measure- 
ments can be performed at HERMES and JLab. 



• Estimations show that agreement with the experi- 
mental data can be improved by means of inclu- 
sion of additional mechanisms which were not in- 
cluded in the model we have presented. 



As a last remark concerning the description of data for 
two-hadron system attenuation in other models. 

There are at least two theoretical works which attempt 
to describe data for two-hadron system NA: first of them 
based on BUU transport model |16], and second based 
on the so called energy loss model 1 17]. 
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APPENDIX 1 

In this appendix we will discuss considerable differ- 
ence between our model and the experimental data at 
the point Z2=0.09 for krypton nucleus, and will try to 
understand the cause of this discrepancy. Our model, 
as mentioned above, takes into account pions only. But 
in heavy nuclei many slow protons are produced, in ad- 
dition to pions (in this discussion we do not distiquish 
kaons from pions). We want to show that by including 
in consideration these protons one can improve agree- 
ment with data. Let us present the two-hadron system 
NA ratio in the form: 



where a is the part of the events which contain pairs 
consisting from a fast pion and a slow proton, R\j and 
R\f are two-hadron system NA ratios in case when a 
pair consist of two pions and pion-proton, respectively. 
It is convenient to introduce a parameter j3 = Rlf /Rm- 
Then Rj^ can be rewritten in the form: 



r>2h 

R M 



(1 + a(0 



(8) 



r>2h 

R M 



(1 



aRjf 



(7) 



Parameter (3 can be defined from the data on single 
hadron attenuation, if one can assume that the fast 
and the slow hadrons in the event are produced inde- 
pendently (correlation between them can be neglected). 
Then (3 w R p M lR\ { , where R P M and B? M are single 
hadron NA ratios for proton and pion respectively In 
addition to the z 2 = 0.09 point we make estimates 
for the next point = 0.15 also, because contribu- 
tion of protons in this point can still be considerable. 
Using available data for krypton |18], we have extrap- 
olated it in the point z = 0.09(0.15) and obtained 
i?£f w 1.5(1.425) and R* M rj 0.85(0.85), which gave 
(3 « 1.765(1.676). Then we determined the parame- 
ter a from eq. (8). For z 2 = 0.09(0.15) the model 
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gives w 0.89(0.87), while the experimental value 
is R$ « 1.2(0.92) which gives a « 0.47(0.087). The re- 
sults of calculations with the model improved by means 
of inclusion of the protons, as was discussed above, are 
shown in Fig. 5a for the krypton case. Three theoreti- 
cal curves are shown: the solid curve corresponds to the 
TSM with CFT1; dashed curve corresponds to the ITSM 
with CFT2; dotted curve corresponds to the ITSM with 
CFT1. Such a value of a for the first point seems to be 
too big. But one can note that still, this correction works 
in the right direction and improves agreement of theory 
with the experimental data even if real value of param- 
eter a is smaller than that following from our estimate. 
In Fig. 5b the dependence on value of parameter a at 
Z2 — 0.09 is shown. As an example we take the case of 
TSM with CFT1 calculated for three values of a equal 
to 0.47 (upper curve), 0.30 (middle curve), 0.10 (lower 
curve). Simultaneously we proportionally changed val- 
ues of a in the next experimental point also. The es- 
timates with other versions of the model give close re- 
sults. 



APPENDIX 2 

As mentioned in text, the Fermi motion of nucleons 
in nucleus and the presence of nucleon-nucleon correla- 
tions (fluctons) can improve agreement with the expe- 
rimental data at Z2 = 0.44 by means of extending the 
range of integration. 

Let us first consider the influence of the Fermi motion 
on increase of z max . We would like to remind the reader, 
that in case of scattering on rest nucleon z max = 1. The 
total centre of mass energy of the secondary hadrons W, 
and in case of taking into account Fermi momentum we 
mark it as Wf- Then, after averaging over the angle be- 
tween virtual photon and nucleon momenta, it can be 
presented as: 

Wp = M 2 - Q 2 + 2E N v = M 2 -Q 2 + 2Mv ef f , (9) 

where M and En are mass and energy of nucleon re- 
spectively, v e f f is the value of v which give on the rest 
nucleon value of W equal Wf- Then for z max we obtain: 



Now we can estimate the influence of the Fermi 
momentum. If we take the average Fermi mo- 
mentum for middle and heavy nucleus to be equal to 
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0.25 GeV/c d, then: z max =1.035 and NA ratio with 
taking into account that also the integration region for 
single hadrons will be increased (about 10%), R^} at 
Z2 = 0.44 must be multiplied on factor 1.42. 

As an alternative mechanism the nucleon-nucleon 
correlations could be considered in the framework of the 
flucton mechanism. That means, that a virtual photon 
scatters on fluctons with the masses M, 2M, 3M, etc. If 
one takes into account only one- and two-nucleon corre- 
lations, one obtains: 

Z m ax = Wl + tt/| , (11) 



where ctfi is the probability of scattering on the two- 
nucleon flucton. If, for the sake of an estimate, one takes 
a ft = 0.01 then z max = 1.01 and the NA ratio Rf} at 
Z2 = 0.44 must be multiplied by a factor of 1.17. We see 
that both mechanisms give considerable improvement 
for the last z 2 point. 
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